On a weak solution of Einstein equations for expanding dust 
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An expanding spherically symmetric dust cloud is considered in a framework of general relativity. 
Initial conditions leading to a shell-crossing singularity are chosen. The way to construct a weak 
solution for such a case is proposed. Suggested method consists in cutting off the region containing 
the shell- crossing and matching the remaining parts of space-time at a thin shell. Junction conditions 
determine the motion of that thin shell. The singular part of dust stress-energy tensor is nontrivial 
only after the shell-crossing occurs. Before that the solution coincides with Lemaitre - Tolman - 
Bondi one. A toy model representing an underdensed region in Universe is discussed. 
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PACS numbers: 04.20.Dw, 04.20.Ex 

The Lemaitre - Tolman - Bondi (LTB) solution de- 
scribing the evolution of a dust cloud has recently be- 
come quite popular for studying inhomogenities in the 
Universe. One of the reasons behind that is the oppor- 
tunity to explain accelerating expansion of the Universe 
without introducing dark energy. See [l| for a compre- 
hensive review and clarification of the idea (before 2008) 
and for a more recent progress. Another important 
feature of the LTB space-time is a possibility of describ- 
ing voids formation (see [s^] including a historical review 
of voids discovery). And even more cosmological and the- 
oretical applications of that remarkable solution can be 
found in [J|. 

One of the less studied properties of the LTB solution 
is the formation of shell-crossing singularities (SCS) for 
certain initial conditions. The cause for it is the intersec- 
tion of initially different dust layers resulting in diverging 
and even negative density. The employing of frameworks 
other then a co-moving one merely brings the metric ten- 
sor to a regular form [5| but can't remove the singularity 
because it is in fact a physical but not a coordinate ef- 
fect. For that reason the initial conditions leading to a 
SCS are usually avoided even if it seems unfortunate. 

The nature of SCS was investigated by different au- 
thors and the conclusion is that it has a different 
("weak" or "inessential") type from a shell- focusing sin- 
gularity and therefore the solution can be extended be- 
yond the SCS. The first example of such an extension 
was provided in Q for a rather special case of space- 
time. Further works Jsl, [l^ suggest the extension to be 
a weak sokition of Einstein equations. In j7l] | such weak 
solutions are derived treating SCS as a shock wave and 
using Rankine - Hugoniot conditions. Unlike the classi- 
cal solution the weak one is not unique to the future of 
the shell-crossing singularity even for well posed initial 
conditions. There is a weak or extended solution which 
has singular part in stress-energy tensor and there is still 



a classical solution which is a special case of a weak so- 
lution with only regular distributions involved. 

Here we introduce another way to find a weak solution 
employing^ Israel - Darmois - Lichnerowicz junction for- 
malism [l2| . The idea is to cut out the unphysical regions 
with negative density and match the remaining parts of 
the space-time at a thin shell. Of course models with 
thin shells are nothing new in cosmology and were first 
studied in [l^. However the important difference is that 
the thin shell in present work arises from smooth initial 
conditions. 

Israel - Darmois - Lichnerowicz matching procedure 
can be viewed as a consequence of dealing with field equa- 
tions in a framework of tensor distributions [1^ - The 
same is true for relativistic Rankine - Hugoniot equa- 
tions which are identical to O'Brien - Singe conditions 
while written in general form of energy and momentum 
conservation across the junction surface [Tsj . Thus re- 
lying on the matching scheme one can expect to get the 
same results avoiding complicated issues of applying gen- 
eralized functions to nonlinear theory. 

The units with G — c = 1 are used throughout. 



I. JOINING TWO LEMAITRE - TOLMAN - 
BONDI SPACE-TIMES AT A THIN SHELL 

In comoving coordinates the line element for dust LTB 
solution is 



r'\T,R) 
1 + E{R) 



dR^ -r^{T, R)dVi^ 



(1) 



where dil 



de^ 



Each layer of the dust 
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is marked with it's own value of radial coordinate R. 
Function r(T, R) determines the distance to the cen- 
ter of the dust particles with given R. Shell-crossing 
singularity appears when the layers of dust intersect 
each other and so the particles with different R start 
to have the same value of radial function r{T,R). This 
leads to a multi- valued behavior of Misner - Sharp mass 
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m(T, r) = m(i?) = (f^ — ii^) r/2 as a function of r [Tl]. The matching yields 
To get rid of the ambiguity one can cut out the lay- 
ers with multi- valued mass from some point R — Ri (r) r[ / dRi 
till R — R2{t) (see Fig. [T|). Functions Ri{t) and R2{t) 1 + Ei \ dr 
should then be found by matching the remaining interior 
and exterior. 
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FIG. 1: Radial metric function r{T,R) after the formation of 
SCS (the expanding interval [Ri(t), R2{t)] should be cut off 
and replaced with the thin shell), T2 > ri 

Not a boundary surface but a thin shell should be 
used while matching because of the gap between masses 
m(i?i) and m{R2). The equation of this shell is specified 
by R — Ri (t) for the inner part of dust cloud and by 
R = R2{t) for the outer part. The intrinsic coordinates 

are {r, 6, 0} for both. The unit normal has following 
components 
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The first equation has two solutions but only the one 
with opposite signs provides the required mass disconti- 
nuity 



dRi 



y/TTE{ dr ^/TT^h dr 

The other solution gives only the classical LTB space- 
time. 

The full time derivative of ([6]) 



, dRi , , o?i?2 
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together with ([7|) forms a linear system on dRi/dr and 
dR2/dT. The solution of this system is 



dRi^ 
dr 



7-2,1 - ^■1,2 



(9) 



The velocity of the thin shell follows immediately from 
the above equation 



dr 



dRia rl^/l + El +7-2^/1 + E2 



VI + £^1 + Vi + E2 



(10) 



For marginally bound case this coincides exactly with the 
shock velocity derived in (ll| from Rankine - Hugoniot 
equations. 



II. A TOY MODEL 
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where index 1 denotes the values of all functions at 
R = Ri (r) while index 2 does the same for R = R2 (r) . 
The only junction condition is continuity on the shell of 
the hrst fundamental form with the components 



gab = gapela)el^) = 9af3-g^-g^- 
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For interior and exterior parts of the considered dust it 
is given by 
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Let's consider an LTB space-time with a negative cur- 
vature as an example. The solution has the parametric 
form 



m{R) 
E{R) 
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(cosh 77 — 1) , 

(sinh77-r;)-ro(i?). (11) 
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With smooth initial conditions one can always choose 
r(r = 0, i?) = R. With that the solution depends on only 
two arbitrary functions, say the initial energy density 
and the initial parameter -q prohles po{R), 'i]o{R). All the 
other functions can be expressed as follows 
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m(R) = / 4Trpo{R)R'^ dR, 



(12) 



3 



i?(^)-^(cosh%(i?)-l), (13) 



ro(i?) = (sinhr;o(i?) - Vo{R)) ■ (M) 

The square of initial velocity has the form 
2m(i?) _ m{R) 



R 



R 



(cosh 770 (i?) + 1) . 



_ (15) 
Following [3] let's take initial conditions at the time of 
last scattering. It is shown in Q that final state of the 
evolution is less sensitive to the initial density perturba- 
tions then to the initial velocity perturbations. So for the 
sake of simplicity we use a homogeneous initial density 
profile po — const. The shape of 770 (i?) function is chosen 
in a way that allows formation of a SCS, namely 

?7o(i?) -ryoc + (?7c-r/oc)(l + a2i?V^o)e-"^'^'/'^«, (16) 

where rjc and rjao determine the Friedmann-likc behavior 
of the curvature near the center and at infinity conse- 
quently: 
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unique solution si = 82- An arbitrary co-moving surface 
can be chosen as a junction surface so one can set also 
si — S2 = Sscs- But this single root splits into several 
distinct ones as the shell-crossing occurs. For smooth 
initial conditions the above splitting first appears when 







and 







(18) 



simultaneously. 

Instead of dealing with the DAE system ^ and (0 
directly one can use standard methods to solve it's con- 
sequence (jlOp which is merely an ODE to find the mo- 
tion of the thin shell. Ri{t) and R2{t) are then calcu- 
lated from r(r) = r{T, Ri{t)) = r{T,R2{T)). The results 
are displayed at Fig. [51 

Present time Hubble constant and density of the weak 
solution are shown at Fig. [31 Vertical dotted lines corre- 
spond to a positions of the thin shells. Both interior and 
exterior parts of the weak solution coincide with classical 
LTB solution but there is a gap between them. 
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FIG. 2: Numerical solution of junction conditions with r\ao ~ 
0.01, 77c = 0.15, ai = 15, 02 = 20 

Now we should solve the system of matching conditions 
^ and ([7]). For numerical calculations it is convenient 
to work with a new variable s — R/ Rq The first step is to 
find a point {tscs, Sscs) of the globally earliest occurrence 
of the shell-crossing singularity. Before the shell-crossing 
formation the considered matching is trivial and © has a 
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FIG. 3: Hubble constant and density at the present time: 

1 — r?oo = 0.01, 77e = 0.15, ai = 15, aa = 20, 

2 — 7700 = 0.01, 77c = 0.09, ai = 10, 02 = 30 
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III. DISCUSSION 

Stress-energy tensor of the extended solution has the 
form 

T''' = T^L + S^'efy^^^S^{^). (19) 

Here (5s (x) is a Dirac's deha function with support on the 
thin sheU. T^^^^^ is a stress-energy tensor of the dust. S""^ 
is a surface stress-energy tensor related to the extrinsic 
curvature K"''' via the Lanczos equation 

Sab = -^{[Kat]~[K:]gat). (20) 

Before the SCS occurs the components of Sab are all equal 
to zero and the space-time is described by Lemaitre - 
Tolman - Bondi solution. After the SCS appears all di- 
agonal components of Sab become nontrivial and can now 



be expressed in the form of a perfect fluid 

Sab = {cr + p)VaVb - Pgab (21) 

where Va = Vae"^-^ and the 4- velocity of the shell differs 
from the 4-velocity of the dust. Because of the spherical 
symmetry the expressions for surface energy density and 
pressure are simply a = S^ and p = —Sg = —S'lj^. Both 
of them are positive in our toy model. 

Thereby we have constructed a weak solution of Ein- 
stein equations with the same initial conditions as classi- 
cal LTB solution and with positive energy density every- 
where including the thin shell. So the main conclusion 
for that paper is that initial conditions leading to shell 
crossings should not be forfeited just because of the sin- 
gularity as it can be replaced with a thin shell. Such 
initial conditions should be considered equally with any 
others. 
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